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Abstract 

A  model-based  electrochemical  impedance  spectroscopy  (EIS)  approach  that  combines  an  equivalent  electrical  circuit  (EEC)  method  and  a 
mathematical  model  derived  from  the  reaction  kinetics  is  proposed  to  investigate  the  simultaneous  oxygen-reduction  reaction  (ORR)  and  methanol- 
oxidation  reaction  (MOR)  at  the  cathode  of  a  DMFC.  Good  agreements  between  the  calculated  results  and  the  experimental  data  validated  the 
proposed  method.  Detailed  kinetic  parameters  and  state  variables  of  the  cathode  were  conveniently  extracted  and  the  concerned  reaction  processes 
were  further  analyzed,  which  demonstrated  the  comprehensive  applicability  of  this  method.  The  results  showed  a  significant  poisoning  effect  on  the 
ORR  by  the  presence  of  methanol  at  the  cathode.  The  results  also  indicated  that  whether  the  methanol  permeated  from  the  anode  can  be  completely 
oxidized  by  electrochemical  reaction  at  the  DMFC  cathode  depends  on  the  electrode  design  and  operating  conditions,  even  at  high  potentials. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Methanol  crossover,  the  permeation  of  methanol  through  the 
ionomer  membrane  from  the  anode  to  the  cathode,  results  in 
simultaneous  cathodic  oxygen-reduction  reaction  (ORR)  and 
anodic  methanol-oxidation  reaction  (MOR)  at  the  cathode  of 
a  direct  methanol  fuel  cell  (DMFC).  The  simultaneous  MOR 
and  ORR  lead  to  not  only  a  fuel  loss,  but  also  a  decrease  in  the 
cathode  potential,  which  is  one  of  major  problems  for  DMFC 
performance  [1-3].  To  suppress  the  effect  of  the  permeated 
methanol,  it  is  essential  to  gain  better  understanding  of  the  mech¬ 
anism  of  how  the  MOR  interacts  with  the  ORR  and  how  this 
interaction  lowers  the  cathode  potential.  Recently,  to  understand 
the  mechanism  of  the  simultaneous  MOR  and  ORR,  extensive 
investigations  have  been  reported  by  using  traditional  electro¬ 
chemical  and  spectroscopic  methods  [1-8].  Despite  a  significant 
progress  made,  the  kinetics  of  simultaneous  MOR  and  ORR  at 
the  DMFC  cathode  is  still  not  well  understood,  mainly  due  to  the 
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complicated  nature  of  the  reactions  involving  multi-step  kinet¬ 
ics,  and  hence  further  studies  are  needed. 

Electrochemical  impedance  spectroscopy  (EIS)  is  a  powerful 
tool  for  elucidating  reaction  mechanisms  involving  multi-step 
kinetics,  because  it  has  the  advantage  of  separating  different 
rate  processes  in  the  frequency  domain,  and  therefore  provides 
a  better  insight  on  the  interfacial  processes  and  species  in  electro¬ 
chemical  systems.  Up  to  now,  however,  only  little  attention  has 
been  paid  to  the  application  of  EIS  to  the  investigation  of  the 
DMFC  cathode  behavior  [9-13].  Muller  and  Urban  [9]  firstly 
introduced  a  method  of  obtaining  DMFC  cathode  impedance 
spectra  based  on  recording  individually  the  impedance  data  for 
the  entire  DMFC  and  the  anode  only,  followed  by  the  calculation 
of  the  DMFC  cathode  impedance.  Using  the  same  procedure, 
Piela  et  al.  [10]  measured  and  gave  some  further  interpretation  of 
the  cathode  impedance  spectra.  Also,  aiming  to  examine  trans¬ 
port  and  reaction  kinetics,  some  reference  electrode  methods 
were  adopted  to  measure  the  cathode  impedance  spectra  individ¬ 
ually  [11-13].  These  pioneering  EIS  investigations  are  greatly 
helpful  in  understanding  the  cathode  behavior.  However,  partly 
because  of  the  complexity  of  the  impedance  data,  all  these  inves¬ 
tigations  used  the  equivalent  electrical  circuit  (EEC)  method,  in 
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which  the  obtained  parameters  may  not  have  physical  equiva¬ 
lents,  and  therefore,  it  is  difficult  to  extract  detailed  information 
about  the  electrode  reaction  mechanism  [14].  In  addition,  for  a 
given  impedance  spectra,  there  may  be  several  EECs  that  can 
all  provide  a  good  fitting,  making  it  difficult  to  select  the  proper 
EEC  model  for  the  studied  electrode  system. 

In  this  work,  we  propose  a  model-based  EIS  approach,  which 
combines  the  EEC  method  and  a  mathematical  model  estab¬ 
lished  from  the  detailed  reaction  mechanism,  and  apply  this 
approach  to  investigate  the  simultaneous  MOR  and  ORR  at  the 
DMFC  cathode.  This  method  can  avoid  the  problems  of  the  only 
EEC  method,  and  more  importantly,  provide  very  detailed  reac¬ 
tion  information,  e.g.  the  kinetic  parameters  and  state  variables. 
To  the  best  of  our  knowledge,  this  is  the  first  attempt  to  study 
the  DMFC  cathode  by  a  model-based  EIS  method.  We  started 
with  a  mathematical  derivation  of  the  Faradaic  impedance, 
based  on  a  simplified  reaction  mechanism  for  the  simultane¬ 
ous  MOR  and  ORR,  and  established  a  relationship  between  the 
Faradaic  impedance  and  a  general  EEC.  Then,  the  experimental 
impedance  data  were  measured  by  a  half-cell  method.  Finally, 
from  the  established  relationship  and  the  measured  experimen¬ 
tal  data,  the  specific  kinetic  parameters  and  state  variables  for 
the  simultaneous  MOR  and  ORR  were  extracted  and  the  results 
were  analyzed  based  on  the  extracted  parameters. 

2.  Theoretical 


where  0/  and  0/  are  respective  adsorbate  coverage  and  its  change 
rate  with  time. 

Apparently,  to  determine  Zp  as  a  function  of  co  at  a  fixed 
potential  E ,  the  relationship  between  7p,  0 ;,  (9Z- ,  and  E  has  to  be 
identified  to  extract  the  parameter  values  for  impedance  calcula¬ 
tion.  To  this  end,  a  simple  reaction  mechanism  for  simultaneous 
MOR  and  ORR  is  employed: 

Pt  +  CH3OH  Apt  -  COads  +  4H+  +  4e“  (10) 

Pt-C0ads+H20Apt  +  C02  +  2H+  +  2e-  (11) 

Pt  +  02  +  H++e-Apt-02Hads  (12) 

Pt  -  02Hads  +  3H+  +  3e  A4Pt  +  2H20  ( 1 3) 

k- 4 

The  electro-catalytic  oxidation  of  methanol  involves  two 
series  steps  [19,20].  The  first  step  is  the  adsorption  and  dehydra¬ 
tion  of  methanol  to  form  the  absorbed  CO  species.  The  second 
step  represents  the  irreversible  oxidation  of  the  CO  species,  lead¬ 
ing  to  CO2  evolution.  Oxygen  reduction  on  Pt  under  acidic  media 
is  usually  characterized  by  the  two-step  Damjanovic  mechanism 
[21,22].  The  first  step  is  an  electrochemical  step  corresponding 
to  the  protonation  of  the  O2,  which  is  the  rate-determining  step. 
Then  the  formed  adsorption  species,  O2H,  is  further  reduced  to 
water.  The  reaction  rates  v;  (i  =  1,  2,  3,  4)  corresponding  to  the 
Eqs.  (10)— (13)  can  be  expressed  as: 


As  discussed  below,  simultaneous  MOR  and  ORR  that  occur 
at  the  cathode  of  a  DMFC  involve  two  adsorbed  species.  Based 
on  the  state  variable  theory  [15,16],  the  Faradaic  impedance  (Zp) 
for  such  an  electrode  process  can  be  expressed  as  [17,18]: 

1 

ZF  “  (l/Rt)  +  (A  +  j a>B)/(D  -  ft)2  +ia>T)  ' 1 ' 

where  co  is  the  radial  frequency, 


1  3/p 

Rt  ~  ~dE 


(2) 


with  7p  and  E  representing  the  Faradaic  current  density  and  the 
electrode  potential,  respectively, 


A  =  m\J\2n2  +  ^2/21^1  —  m  1/22^1  —  ^2^11^2 

(3) 

B  =  m\ni  -f-  m2^2 

(4) 

D  =  J11J22 

-  712/21 

(5) 

T  =  -Jn- 

J22 

(6) 

with 

37p 

mi  =  — 

30/ 

i  =  1,2 

(7) 

30/ 

Hi  =  —  i 

3  E 

=  1,2 

(8) 

J  dbi 

ik  30* 

i,k=  1,2 

(9) 

vi  =  k\c\{\  -  0i  -  02)  exp  (  — 


(  E 

Vl  =  *201  exp  (  — 


v3  =  &3<a(1  -  01  -  02)  exp  (  - 
E\ 


V4  =  k+402  exp 


b\  ) 


03, 

)  -  0-4(1  -  01 


02)  exp 


(14) 

(15) 

(16) 
(17) 


where  k[  (i=  1,  2,  3,  +4,  —4)  are  rate  constants,  bi  (i=  1,  2,  3, 
4)  the  Tafel  slopes,  0i  and  02  the  coverage  of  CO  and  O2H,  ci 
and  C2  the  concentrations  of  methanol  and  oxygen.  Based  on 
the  proposed  reaction  mechanism  (Eqs.  (10)— (13)),  the  Faradaic 
current  density  (7p)  can  be  written  as: 


7p  =  F(4v\  +  2v2  —  V3  —  3v4_) 


(18) 


where  F  is  the  Faradaic  constant.  According  to  mass  balance,  the 
change  rates  of  absorbate  coverage  (0;)  is  the  net  result  of  pro¬ 
duction  and  consumption  of  the  absorbates,  and  can  be  expressed 
as: 

d0i  F 

0\  =  —  =  —  (vi  -  v2)  (19) 

At  q  1 

A62  F 

02=-f  =  ~(V3  -  v4)  (20) 

At  q2 

where  q\  (420  |xC  cnr2  [14])  and  q2  (assumed  420  |xCcm-2) 
are  the  charges  required  for  complete  CO  and  O2H  adsorption 
on  unit  Pt  surface. 
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By  substituting  7p,  Oi  and  0;  into  Eqs.  (2)-(9)  with  the  expres¬ 
sions  in  Eqs.  (14)-(20),  the  relationships  between  Rt,  A,  B ,  D, 
T  and  the  kinetic  and  state  parameters  (kt,  hi ,  Oi  and  cf)  can  be 
established. 

According  to  Wu  et  al.  [18],  the  Faradaic  impedance  for 
an  electrode  process  with  two  adsorbed  species  can  also  be 
expressed  by  a  general  EEC,  shown  in  Fig.  la,  and  thus,  the 
mathematical  expression  of  the  Faradaic  impedance  can  be 
equivalently  written  as 


Zp  —  Rt 


£ 


Ri 

1  +  j  coRiCi 


7  =  1,2 


(21) 


where  R[  is  the  resistance  elements,  Q  the  capacitance  ele¬ 
ments.  Comparing  the  Eq.  (1)  with  the  Eq.  (21),  the  relationship 
between  the  electrical  elements  in  the  EEC  and  the  electrochem¬ 
ical  parameters  in  the  Faradaic  impedance  expression  Eq.  (1)  can 
be  established  as: 


A=l(  *1  +  *2  N 

R2,  \RlR2ClC2 ) 

1 

D  = - RtA 

R1R2C1C2 


(22) 

(23) 

(24) 


R\C\  +  R2C2 

R1R2C1C2 


-  RtB 


(25) 


From  Eqs.  (22)-(25),  the  electrochemical  parameters  A,  B ,  D  and 
T  can  be  calculated  from  the  values  of  the  electrical  elements 
Ri ,  Ri  and  Q,  which  can  be  obtained  by  fitting  experimental 
data  with  the  proposed  EEC.  Also,  as  was  mentioned  above, 
A,  B ,  D  and  T  as  well  as  Rt  are  related  to  the  kinetic  and  state 
parameters  (£;,  Oi  and  cf)  of  the  simultaneous  MOR  and  ORR 


kinetics.  Therefore,  based  on  Eqs.  (2)-(6),  the  kinetic  and  state 
parameters  (&;,  bu  Oi  and  q)  can  be  obtained  from  the  calculated 
A,  B ,  D ,  T  and  Rt  by  a  simplex  nonlinear  regression  procedure. 

3.  Experimental 

The  experimental  impedance  spectra  of  the  DMFC  cathode 
were  measured  by  a  half-cell  method,  which  was  proved  reliable 
for  EIS  measurement  [23-25].  In  this  half  cell,  a  half  membrane 
electrode  assembly  (MEA),  consisting  of  a  Nation  membrane 
with  a  gas  diffusion  electrode  (GDE)  on  one  side,  was  held 
vertically.  The  GDE,  consisting  of  a  catalyst  layer  and  a  gas 
diffusion  layer,  was  in  contact  with  a  flow  field  plate,  through 
which  O2  was  supplied  and  electric  current  was  extracted.  The 
membrane  was  exposed  to  the  electrolyte  solution,  in  which  a 
platinum  mesh  was  immersed  to  serve  as  a  counter  electrode. 
A  saturated  calomel  electrode  (SCE)  was  positioned  approxi¬ 
mately  4.0  mm  away  from  the  membrane  via  a  Luggin  capillary 
to  serve  as  a  reference  electrode.  The  half-MEA  was  obtained  by 
hot-pressing  a  Nation  112  membrane  and  a  1  cm2  in-house  GDE 
with  Pt  loading  0. 1  mg  cm-2.  The  in-house  fabricated  GDE  was 
prepared  by  spraying  a  catalyst  ink  consisting  of  60  wt.%  Pt/C 
catalyst,  5  wt.%  Nation  solution  and  ethanol  onto  a  wet-proofed 
carbon  cloth  substrate. 

All  the  electrochemical  measurements  were  conducted  in 
deaerated  0.5  M  H2SO4  +  O.5M  CH3OH  solution  at  ambient 
temperature  (20  °C),  using  an  Autolab  PGSTAT  30  with  a  fre¬ 
quency  response  analyzer  module,  FRA2.  A  high  flow  rate 
of  unhumidified  O2  (100  seem)  at  atmospheric  pressure  was 
adopted  to  eliminate  the  mass-transport  limitations  of  oxy¬ 
gen  and  keep  the  ORR  almost  under  pure  kinetic  control.  The 
impedance  spectra  were  recorded  10  points/decade  in  the  poten- 
tiostatic  mode  by  sweeping  frequencies  from  50k  to  0.02  Hz 
range  at  an  amplitude  of  10  mV. 


(a) 


(b)  C0  ZCPE 


*0 


Fig.  1.  Equivalent  electrical  circuit  for  the  evaluation  of  impedance  responses 
of  the  DMFC  cathode:  (a)  Faradic  impedance  for  an  electrode  reaction  with 
two  absorption  intermediates,  and  (b)  a  complete  circuit  model  for  the  DMFC 
cathode. 


4.  Results  and  discussion 

Fig.  2  shows  the  experimental  EIS  patterns  of  the  electrode  at 
various  potentials.  At  high  frequencies,  the  impedance  showed  a 
line  with  a  slope  close  to  45°.  This  line,  apparently  independent 
of  the  electrode  potential,  could  be  ascribed  to  the  ionic  ohmic 
drop  in  the  membrane  or  the  electronic  contact  between  different 
layers  [21,26].  At  medium  frequencies,  a  large  capacitive  arc 
was  observed,  which  was  associated  with  the  charge  transfer 
processes  [27]  and  depressed  due  to  a  high  roughness  of  the 
electrode.  At  low  frequencies,  there  was  an  unobvious  pseudo 
inductive  loop.  Such  pseudo  inductive  patterns,  which  are  typical 
features  of  systems  with  adsorbed  intermediates  [28],  are  critical 
to  proposing  the  reaction  kinetics  and  the  fitting  EEC  for  the 
DMFC  cathode. 

Considering  the  measured  impedance  patterns,  the  reaction 
mechanism,  and  the  structural  features  of  the  half  cell,  the  exper¬ 
imental  data  in  Fig.  2  were  fitted  by  a  proposed  EEC  shown  in 
Fig.  lb,  where  Ro  is  the  membrane  or  contact  resistance,  Co 
the  membrane  or  contact  capacitance,  Zcpe  the  constant  phase 
element  (CPE),  and  Zp  the  Faradaic  impedance  as  presented  in 
Fig.  la.  In  this  EEC,  the  parallel  combination  of  Rq  and  Cq  was 
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Table  1 

Fitted  values  of  the  elements  in  the  equivalent  electrical  circuit  illustrated  in  Fig.  1 


£(V) 

Y  (F  m-2) 

a 

Ro  (£2  m2) 

C0  (Fm-2) 

Rt  (£2m2) 

Ri  (£2  m2) 

Cl  (Fm-2) 

R2  (£2  m2) 

C2  (Fm2) 

0.75 

0.1119 

0.7997 

0.00156 

0.01148 

0.10480 

0.16080 

0.03190 

-0.01412 

-19.01 

0.71 

0.1031 

0.8275 

0.00156 

0.01408 

0.07910 

0.08720 

0.03611 

-0.01030 

-15.73 

0.68 

0.1101 

0.8010 

0.00155 

0.01267 

0.05008 

0.06852 

0.03874 

-0.007580 

-14.92 

0.65 

0.1187 

0.7910 

0.00159 

0.01533 

0.04246 

0.04712 

0.04637 

-0.006090 

-12.68 

Re(Zv)  /  ohm  cm2 


Fig.  2.  Complex-plane  impedance  patterns  of  the  electrode  at  various  potentials. 
Scattered  points:  experimental  data;  solid  line:  the  fitted  data  by  the  equivalent 
circuit  illustrated  in  Fig.  1. 


used  to  fit  the  ionic  ohmic  drop  in  the  membrane  or  the  elec¬ 
tronic  contact  between  different  layers.  The  CPE,  instead  of  a 
pure  capacitor,  was  used  for  the  electric  double-layer  capaci¬ 
tance  to  account  for  the  electrode  roughness.  The  impedance  of 
a  CPE  can  be  written  as 

Zcpe  =  .  „  (26) 

Y(jco)a 

where  Y  is  the  magnitude  and  a  the  dimensionless  parameter 
which  has  a  value  from  —1  to  1.  The  fitted  results  using  the 
above  EEC  are  listed  in  Table  1  and  the  simulated  EIS  patterns 
with  these  fitted  results  are  shown  in  Fig.  2,  in  which  a  good 
agreement  with  the  experimental  data  was  observed.  Based  on 
the  established  correlations  (Eqs.  (22)-(25)),  the  electrochem¬ 
ical  parameters  A,  B ,  D  and  T  in  the  mathematical  expression 
of  Faradaic  impedance  (Eq.  (1))  at  various  potentials  were  then 
calculated  from  the  fitted  values  of  the  electrical  elements  and 
are  given  in  Table  2. 

After  obtaining  the  electrochemical  parameters  in  the  Farad¬ 
aic  impedance  expression,  the  specific  kinetic  parameters  and 
state  variables  for  the  simultaneous  MOR  and  ORR  mechanism 


Table  2 

Calculated  values  of  the  electrochemical  parameters  in  Eq.  (1) 


£(V) 

A  (Q  1  m  2  s  2) 

5(f2-1m-2  s— 2) 

D  (s-2) 

r(s-‘) 

0.75 

-9.70  x  103 

-2.85  x  103 

1.74  x  103 

497.3 

0.71 

-2.17  x  104 

-3.97  x  103 

3.48  x  103 

606.5 

0.68 

-8.09  x  104 

-1.03  x  104 

7.38  x  103 

899.7 

0.65 

-1.35  x  105 

-1.19  x  104 

1.17  x  104 

976.7 

Table  3 

Fitted  kinetic  parameters  for  the  simultaneous  MOR  and  ORR. 


Parameters 

Values 

k\  (mols-1  m-2) 

2.86  x  10-8 

k,2  (mols-1  m-2) 

6.48  x  10-11 

k3  (mols-1  m-2) 

0.28 

k+ 4  (mols-1  m-2) 

88.70 

k-4  (mols-1  m-2) 

3.96  x  10-8 

bi  (V  dec-1) 

0.058 

b2  (V  dec-1) 

0.069 

b3  (V  dec-1) 

0.053 

MV  dec-1) 

0.062 

were  extracted  by  a  nonlinear  fitting  procedure  based  on  the 
established  correlations  (Eqs.  (2)-(6)).  Table  3  summarizes 
the  fitted  kinetic  parameters  for  the  simultaneous  MOR  and 
ORR.  For  MOR,  ki  is  almost  two  orders  of  magnitude  lower 
than  k\ ,  indicating  oxidation  of  the  absorbed  CO  was  the  rate¬ 
determining  step  whereas  adsorption  of  the  methanol  molecule 
to  form  the  CO  species  was  relatively  faster.  For  ORR,  k 3  is 
significantly  lower  than  k+ 4,  indicating  the  oxygen  reduction  to 
the  absorbed  O2H  was  the  rate-determining  step,  which  was  in 
agreement  with  the  previous  assumption  for  the  ORR. 

The  coverage  of  the  absorbed  species  at  a  DMFC  cathode, 
which  is  often  obtained  by  complicated  spectroscopic  appara¬ 
tuses,  could  also  be  obtained  by  our  model-based  EIS  method. 
Fig.  3  presents  the  calculated  values  of  CO  and  O2H  coverage,  as 
well  as  the  free  surface,  as  a  function  of  the  electrode  potential. 
With  the  increase  in  potential,  0\  decreased  whereas  62  increased 
steadily.  From  0.65  to  0.75  V,  0\  decreased  from  about  0.7  to 
0.1  and  62  increased  from  0.3  to  0.8.  The  continuous  declining 
of  6 1  could  be  interpreted  by  the  electro-oxidative  removal  of 


Fig.  3.  Calculated  CO  and  O2H  coverage  (6fi  and  62)  and  free  surface 
(1  —  6\  —  62)  as  a  function  of  electrode  potential. 
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CO  species  at  higher  potentials.  The  increase  of  62  was  prob¬ 
ably  due  to  the  competitive  adsorption  between  O2H  and  CO. 
The  significant  removal  of  CO  at  higher  potentials  led  to  more 
active  sites  available  for  O2H  adsorption,  and  thus,  O2  increased 
with  potential.  Noteworthy  was  the  large  CO  coverage,  espe¬ 
cially  at  lower  potentials,  which  was  associated  with  the  strong 
adsorption  strength  of  Pt  with  CO  and  indicated  that  the  pres¬ 
ence  of  methanol  significantly  reduced  the  active  sites  for  the 
ORR  at  the  DMFC  cathode.  The  free  active  surface,  also  shown 
in  Fig.  3,  was  quite  small,  which  was  consistent  with  the  above 
analysis. 

In  order  to  show  the  methanol  effect  on  the  ORR  more 
straightforwardly,  the  polarization  curves  were  calculated  and 
are  presented  in  Fig.  4.  Clearly,  the  polarization  curve  for  simul¬ 
taneous  MOR  and  ORR  showed  significant  lower  performance 
than  that  for  the  pure  ORR,  which  was  calculated  by  assum¬ 
ing  zero  methanol  concentration  and  CO  coverage,  especially  at 
lower  potentials.  This  result,  consistent  with  the  coverage  varia¬ 
tion,  indicates  that  the  presence  of  methanol  greatly  affected  the 
ORR.  The  total  effect  of  methanol  on  the  ORR  could  be  divided 
into  two  parts:  the  parasitic  current  effect,  which  assumed  MOR 
and  ORR  independent  of  each  other,  and  the  poisoning  effect, 
which  was  mainly  due  to  adsorption  of  the  CO  species.  To  fur¬ 
ther  clarify  the  methanol  effect  on  ORR,  the  parasitic  current 
effect  and  the  total  effect  were  respectively  obtained  from  the 
partial  MOR  current  density  and  the  difference  of  the  current 
densities  between  pure  ORR  and  simultaneous  MOR  and  ORR 
at  the  same  potential.  The  ratio  of  parasitic  current  effect  and 
total  effect  was  plotted  in  the  inset  of  Fig.  4.  The  ratio  of  para¬ 
sitic  current  effect  and  total  effect  was  below  0.1  at  0.65  V  and 
then  increased  to  over  0.5  at  0.75  V,  revealing  that  the  poisoning 
effect  was  dominant  at  low  potentials  while  the  parasitic  current 
effect  was  more  important  at  high  potentials,  which  was  con¬ 
sistent  with  the  results  in  Fig.  3.  These  results  indicate  that  the 
interaction  between  ORR  and  permeated  methanol  could  not  be 
explained  by  a  classical  mixed  potential  theory,  especially  at  low 
potentials. 


Current  density  /  mA  cm*2 

Fig.  4.  Calculated  polarization  curves  for  pure  ORR  and  simultaneous  MOR  and 
ORR  at  various  potentials.  Inset:  the  ratio  of  parasitic  current  and  total  effects 
of  permeated  methanol  on  ORR. 


Methanol  crossover  through  the  membrane  occurs  mainly 
due  to  the  driving  forces  of  concentration  gradients  and  electro¬ 
osmosis,  which,  at  steady  state,  can  be  described  by 


/p  3c  1 

k— — -  =0 
F  dx 


(27) 


where  Dm  is  the  methanol  diffusivity  in  Nation  membrane,  k 
the  coefficient  associated  with  electro-osmotic  drag.  This  equa¬ 
tion  was  solved  to  give  the  methanol  permeation  rates  due  to 
the  concentration  diffusion  and  the  electro-osmotic  drag.  Dur¬ 
ing  the  calculation,  the  values  of  1.37  x  10“ 10  m2  s-1  and  0.045 
were  used  for  Dm  and  A,  which  were  estimated  from  the  limiting 
current  density  for  methanol  crossover  through  the  Nation  mem¬ 
brane,  as  described  by  Ren  et  al.  [7].  Fig.  5  shows  the  calculated 
methanol  fluxes  resulted,  respectively,  from  the  concentration 
diffusion  and  the  electro-osmotic  drag  at  various  potentials.  For 
the  case  of  the  electro-osmotic  drag,  the  methanol  flux  decreased 
with  increasing  potential,  due  to  lower  cathodic  current  densities 
at  higher  potentials.  For  the  case  of  the  concentration  diffusion, 
on  the  contrary,  a  slight  increase  of  methanol  flux  with  potential 
was  observed,  which  could  be  explained  in  terms  of  the  methanol 
concentrations  at  both  sides  of  the  half  MEA.  Apparently,  the 
methanol  concentration  at  the  electrode  side  was  reduced  due 
to  larger  MOR  rates,  while  the  methanol  concentration  at  the 
other  side  of  the  half  MEA  was  almost  maintained  constant,  at 
higher  potentials.  Hence,  the  methanol  flux  due  to  concentra¬ 
tion  gradient  increased  with  potential.  The  increase  of  methanol 
flux  due  to  concentration  gradient  with  potential,  however,  was 
slight  due  to  the  small  variation  of  the  methanol  concentration  at 
the  electrode  side.  It  was  interesting  that,  even  at  low  potentials 
when  cathodic  current  density  was  relatively  larger,  the  methanol 
flux  due  to  concentration  diffusion  was  still  significantly  higher 
than  that  due  to  electro-osmotic  drag,  which  revealed  that  the 
methanol  flux  due  to  the  concentration  diffusion  was  the  domi¬ 
nant  way  for  methanol  crossover  through  the  Nation  membrane. 
It  should  be  noticed  that  this  was  just  the  observation  in  our  case, 
and  methanol  fluxes  due  to  concentration  diffusion  and  electro- 
osmotic  drag  were  highly  dependent  on  the  electrode  design  and 
operation  condition. 


Fig.  5.  Calculated  methanol  fluxes  resulted  from  the  concentration  diffusion  and 
electro-osmotic  drag  as  a  function  of  electrode  potential. 
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Fig.  6.  The  rate  of  methanol  permeation  through  the  membrane  and  the  rate  of 
methanol  oxidation  at  the  cathode  as  a  function  of  electrode  potential.  Inset:  the 
ratio  of  methanol  oxidation  and  permeation  rates. 


The  rate  of  methanol  permeation  through  the  membrane  and 
the  rate  of  methanol  oxidation  at  the  cathode  as  a  function  of 
the  electrode  potential  were  also  calculated  and  are  compared  in 
Fig.  6.  The  rate  of  methanol  permeation  through  the  membrane 
decreased  with  increasing  potential,  which  was  the  net  result  of 
both  concentration  diffusion  and  electro-osmotic  drag  as  shown 
in  Fig.  5.  The  rate  of  methanol  oxidation,  which  is  difficult  to  be 
measured  by  other  electrochemical  methods  but  can  be  easily 
obtained  by  this  model-based  EIS  method,  increased  steadily 
with  potential,  due  to  faster  MOR  kinetics  at  higher  potentials. 
The  inset  of  Fig.  6  presents  the  ratio  of  the  rate  of  methanol 
oxidation  to  that  of  methanol  permeation.  It  is  interesting  to  find 
that,  at  all  potentials,  the  electrochemical  oxidation  rate  was  only 
40-70%  of  the  permeation  rate,  indicating  that  the  permeated 
methanol  could  not  be  completely  oxidized  by  electro-catalytic 
reaction  at  the  DMFC  cathode,  even  at  high  potentials.  The  non- 
electrochemically  oxidized  methanol  may  remain  in  its  origin 
state,  which  was  proved  by  some  spectroscopic  observations 
[29],  or  be  oxidized  by  a  chemical  reaction,  which  was  observed 
by  Vielstich  and  coworkers  [30,8].  No  matter  what  form  the 
nonelectrochemically  oxidized  methanol  was,  the  assumption 


Fig .  7 .  Simulated  EIS  patterns  based  on  the  extracted  kinetic  parameters .  Hollow 
points:  experimental  data,  solid  points:  simulated  plots. 


that  all  permeated  methanol  could  be  completely  oxidized  by 
electro-catalytic  reaction  might  lead  to  some  errors  in  explaining 
experimental  phenomena  or  mathematical  modeling. 

Finally,  to  verify  the  reliability  of  the  model-based  EIS 
approach  and  the  accuracy  of  the  calculated  results,  the  EIS 
patterns  calculated  from  the  extracted  parameters  were  com¬ 
pared  with  experimental  results  and  are  shown  in  Fig.  7.  Clearly, 
agreements  between  the  experimental  and  simulated  results  were 
observed,  demonstrating  the  reliability  of  the  method  by  the 
combining  use  of  EEC  method  and  a  mathematical  modeling. 

5.  Conclusions 

The  simultaneous  MOR  and  ORR  behavior  at  a  DMFC  cath¬ 
ode,  resulted  from  methanol  crossover,  was  investigated  by  a 
model-based  EIS  method  that  couples  an  EEC  and  a  mathe¬ 
matical  model  derived  from  the  detailed  reaction  kinetics.  This 
method  combines  the  advantages  of  the  quantitative  fitting  of 
the  experimental  data  in  the  EEC  method  and  the  exact  physical 
equivalent  for  every  parameter  in  the  mathematical  modeling, 
and  can  be  easily  implemented  to  extract  the  reaction  details  of 
the  studied  electrochemical  systems.  As  far  as  we  know,  this 
is  the  first  attempt  to  study  the  simultaneous  MOR  and  ORR 
at  a  DMFC  cathode  by  a  model-based  EIS  method.  Agreements 
between  experimental  and  theoretical  EIS  data  demonstrated  the 
reliability  of  the  model-based  EIS  method.  Most  of  the  detailed 
reaction  information,  e.g.  kinetic  parameters  and  state  variables, 
about  the  DMFC  cathode  could  be  conveniently  extracted  out  by 
this  model-based  EIS  approach,  which  is  very  helpful  to  under¬ 
stand  the  cathode  behavior.  From  the  obtained  parameters,  some 
significant  facts  related  to  the  DMFC  cathode  were  revealed.  The 
presence  of  methanol  at  the  cathode  could  lead  to  a  significant 
poisoning  effect  on  the  oxygen  reduction  reaction  due  to  a  high 
CO  coverage,  especially  when  the  cathode  operates  at  low  poten¬ 
tials.  In  modeling  the  DMFC  performance,  it  has  usually  been 
assumed  that  all  permeated  methanol  can  be  completely  oxidized 
by  electrochemical  reaction  at  the  DMFC  cathode.  However, 
this  work  reveals  that  this  assumption  is  questionable.  Whether 
the  permeated  methanol  can  be  completely  oxidized  by  electro¬ 
chemical  reaction  depends  highly  on  the  electrode  design  and 
operating  conditions. 
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